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STATIONARY DISTRIBUTIONS OF MULTI-TYPE TOTALLY 
ASYMMETRIC EXCLUSION PROCESSES 

By Pablo A. Ferrari and James B. Martin 
Universidade de Sao Paulo and University of Oxford 

We consider totally asymmetric simple exclusion processes with 
n types of particle and holes (n-TASEPs) on Z and on the cycle 
Zjv. Angel recently gave an elegant construction of the stationary 
measures for the 2-TASEP, based on a pair of independent product 
measures. We show that Angel's construction can be interpreted in 
terms of the operation of a discrete-time M/M/l queueing server; the 
two product measures correspond to the arrival and service processes 
of the queue. We extend this construction to represent the stationary 
measures of an n-TASEP in terms of a system of queues in tandem. 
The proof of stationarity involves a system of n 1-TASEPs, whose 
evolutions are coupled but whose distributions at any fixed time are 
independent. Using the queueing representation, we give quantitative 
results for stationary probabilities of states of the n-TASEP on Zjv, 
and simple proofs of various independence and regeneration proper- 
ties for systems on Z. 

1. Introduction. Consider first the totally asymmetric simple exclusion 
process (TASEP) on Z or on the cycle Zjv- Some sites iGZ (or Z^r) contain 
a particle; at the others we say there is a hole. The dynamics of the system 
are as follows: a bell rings as a Poisson process of rate 1 at each site indepen- 
dently; when the bell at site i rings, if there is a particle at site i and a hole 
at site i — 1, they exchange. Put another way, each particle in the system 
tries to jump to the left at rate 1; the jump succeeds whenever the site to 
its left is unoccupied. See, for example, [9, 10] for a wealth of information 
on the TASEP and its relatives. 

We will denote a configuration of the TASEP by u = (u(j),j € Z or Zjv) 6 
Ui or U[ N \ where U x = {l,oo} z and u[ N) = {l,oo}V We set u(i) = 1 if 



Received February 2005; revised June 2006. 

A MS 2000 subject classifications. Primary 60K35; secondary 82C22, 90B22. 
Key words and phrases. Totally asymmetric simple exclusion process, multi-type pro- 
cess, multiclass queue. 

This is an electronic reprint of the original article published by the 

Institute of Mathematical Statistics in The Annals of Probability. 

2007, Vol. 35, No. 3, 807-832. This reprint differs from the original in pagination 

and typographic detail. 



1 



2 



P. A. FERRARI AND J. B. MARTIN 



there is a particle at i and u(i) = oo if there is a hole at i. (This notation is 
not standard but will be convenient later.) 

The stationary measures for the TASEP are as follows. On Z^v, the num- 
ber of particles in the system, p, is conserved by the dynamics. For fixed 
p, the dynamics are given by an irreducible Markov chain with finite state- 
space, whose stationary measure is uniform over all ( ) possible states. Any 
stationary measure for the system is some linear combination of these uni- 
form distributions. 

On Z there is a one-parameter family of translation-invariant extremal 
stationary measures v\, A € [0, 1]. Under v\, each site i is occupied by a par- 
ticle independently and with probability A; the measure v\ is the product 
measure with density A of particles (or the measure of a Bernoulli process 
with rate A). The only other extremal stationary measures are the so-called 
blocking measures: the measure concentrated on the configuration with parti- 
cles at all negative sites and holes at all nonnegative sites, and its translates. 

We now consider a TASEP with n classes of particle or n-TASEP, for 
n > 1. Now a configuration u = (u(j),j € Z or Zjy) of the system is a member 
of U n or , where 

Z4 = ({l,2,...,n}U{oo}f, z4 Ar ) = ({l,2,...,n}U{oc})^. 

If u(i) = r <n we say that there is a particle of class r at site i, while if 
u(i) = oo we say that there is a hole at i. 

The dynamics are as follows. Bells ring at rate 1 at each site independently. 
When the bell at site i rings, the values at i — 1 and i swap if u(i — 1) > u(i), 
and remain unchanged otherwise. That is, the process jumps from u to the 
configuration u^ l ~ 1 ^ where 

u^ i) {j) = u{j) i£j<£{i-l,i}, 

U (*-M)(i _ i) = min{n(i - 1), «(<)}, 

u^~ 1,l \i) = max{u(i - l),u(i)}. 

Put another way, each particle in the system tries to jump to the left at rate 
1. The jump succeeds when the site to its left contains a hole or contains 
a particle with a higher class (in which case the two particles exchange 
positions). 

(Note that the way we have defined it, the n-type TASEP has n + 1 
possible values for a site: n types of particle and holes. It might be more 
elegant simply to regard holes as particles of type n + 1, and to replace 
the value oo by n + 1 throughout; but the distinguished state oo will be 
convenient later for various reasons.) 

One natural way in which n-type TASEPs arise is from couplings of sev- 
eral 1-type TASEPs. Let rj[ , rff^ , . . . , rj^ , t > be n processes realizing 
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a TASEP, started from initial conditions such that % (j) — Wo (J) — ' " — 
Wo^H) f° r au 3- (I- e -j whenever there is a particle at site j in process m, 
there is also one in the processes m + 1, m + 2, . . . , n.) Suppose one couples 
the processes by using the same processes of bells at each site for all of 

them. Then the ordering r]\ (j) > T]t 2 \j) > • ■ • > Wt^ii) continues to hold 
for all t (this is an instance of the basic coupling of Liggett [8, 9]). Let 

ut(j) = inf{m : ?][ m ^ = 1} [with Ut(j) = oo if ?][ m ^ (j) = oo for all m]. Then m 
realizes an n-type TASEP. 

The 2-TASEP has been studied from several perspectives. The existence 
of a translation-invariant stationary measure for the process on Z, with 
densities Ai and A2 of first- and second-class particles (for < Ai < Ai + 
A2 < 1), was proved by Liggett [8] to demonstrate ergodic properties of 
the TASEP. The uniqueness and extremality of this measure were shown 
by Ferrari, Kipnis and Saada [5] and Speer [12] (the only other extremal 
invariant distributions are those concentrated on a single state — "blocking 
measures" as for the 1-TASEP above). Derrida, Janowsky, Lebowitz and 
Speer [2] and Speer [12] construct the measure explicitly (both in finite and 
infinite volume) using a matrix method and show various regeneration and 
asymptotic properties; Ferrari, Fontes and Kohayakawa [4] give probabilistic 
interpretations and proofs of the measure and its properties. 

Recently Angel [1] gave an elegant construction of this stationary measure 
based on two independent product measures with densities Ai and Ai + 
A2 (and an analogous construction for the case of Z^v); the proof involves 
providing bijections between certain families of binary trees and pairs of 
binary sequences. 

We will show that Angel's construction can be rewritten in terms of the 
operation of a queueing server (namely, an M/M/l queue in discrete time). 
The two independent product measures correspond to the arrival process 
and the service process of the queue. The stationary measure for the 2- 
TASEP corresponds to the output process of the queue. Sites of the TASEP 
are interpreted as times in the queueing process. If the queue has a departure 
at time i, then one puts a first class particle at site i. If the queue has an 
"unused service" at time i (i.e., a service is available but there is no customer 
present) then one puts a second-class particle at site i. If there is no service 
available at time i, then there is a hole at site i. (Full descriptions are given 
in the next section.) 

We then generalize this result to the n-TASEP with n > 3. First a remark 
about the set of stationary measures. Let Ai, . . . , A n G (0, 1) with J2r=i < 
1. For the process on Z, there is a unique translation-invariant stationary 
measure with density A r of particles of type r, for each r. These stationary 
measures are extremal, and the only other stationary measures are blocking 
measures, concentrated on a single configuration. These facts can be proved 
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almost identically to the proofs given in, for example, [5] and [12] for the case 
n = 2. The case of the n-TASEP on Zat is analogous; for any pi, . . . ,p n 6 
{0,1,..., N} with J2r=iPn < N, there is an extremal stationary measure 
concentrated on configurations with exactly p r particles of type r, for each 
r; these are the only extremal stationary measures. 

We construct a representation of these stationary measures based on a 
system of n — 1 queues in tandem. Take the case of Z, say. Consider a set of 
n independent Bernoulli processes (product measures), such that the ?nth 
has density Ai + • — h A m , for each m. The service process of the mth queue 
is given by the (m + l)st product measure. The first product measure acts as 
the arrival process to the first queue. Thereafter, the arrival process to the 
mth queue, m > 1, is given by the output process from the (m — l)st queue. 
The mth queue is a "priority queue" with m types of customer; at a service 
time, the customer who departs is the one whose class number is lowest out 
of those present. We will show that the distribution of the output process 
of the (n — l)st queue provides the stationary distribution of the n-TASEP 
with densities Ai, . . . , A n . (We describe the operation of the system fully in 
the next section.) 

The proof of stationarity has the following structure. First we define dy- 
namics on a set of n configurations of particles on the line (the "multiline 
process"). On each line, the local transitions are those of a TASEP, and the 
lines are coupled in such a way that the bells on neighboring lines always 
ring at the same time, either at the same site or at neighboring sites. We 
will show that a collection of n independent product measures is stationary 
under these dynamics. (In fact, it will turn out that in this equilibrium, the 
multiline process provides a coupling of n TASEP processes, whose marginal 
distributions at any fixed time are independent.) The second step of the proof 
is to show that, under these dynamics, the final output process derived as 
in the previous paragraph is realizing an n-TASEP. 

In the case n = 2, the two lines of the multiline process correspond to the 
arrival and departure processes of the queue. (We note that, especially in 
the case n = 2, the structure of the proof is similar to that used in recent 
work by Duchi and Schaeffer [3] , where a related two-line process is used to 
study combinatorial aspects of TASEPs with one or two types of particle, 
on Ztv or on a finite strip with boundaries.) 

In Section 2, we give precise definitions and state the main results, first for 
n = 2 (Angel's construction and the queueing interpretation) and then for 
general n. In Section 3, the first step of the proof is carried out: the "multiline 
process" is defined and it is proved that a collection of independent product 
measures is stationary under its dynamics. (For the case of Zjv, the analogous 
stationary distribution is a collection of independent uniform distributions 
over configurations with a fixed number of particles.) The second step of the 
proof described in the previous paragraph is done in Section 4. 
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In Section 5 we use the multiline representation to prove a quantitative 
result about the probabilities of states in the stationary distribution of an 
n-TASEP on the ring Z N . Let p u . . . ,p n e {0, 1, . . . , A} with £™ =1 p r < A 
and consider the stationary distribution concentrated on configurations with 
exactly p r particles of class r for each r. Then the probability of every state 
is an integer multiple of M~ l , where 



(1.1) M 



N\f N \ t A 

Pi) KPl+Pz) "' \Pl +P2 H VPn 



This confirms a conjecture of Angel [1]. (See also the discussion by Maffick, 
Mallick and Rajewsky in [11], where the matrix method is used to investigate 
the case n = 3.) We can then identify those states with probability M _1 (the 
smallest possible). A configuration u is minimal in this way if the following 
holds for each j € 1*n'- if u(j) = oo then u(j + 1) = oo or u(j + 1) = n, while 
if u{j) = m < oo then u(j + 1) > m — 1. 

For example, consider n = 2, = 9, p\ = 3, p2 = 3. Possible "minimal 
states" include: (oo, oo, oo, 2, 2,2, 1, 1, 1) and (oo, 2, oo, oo, 2, 1, 2, 1, 1). 

For the extreme case of N types of particle on Zjy, with p m = 1 for all 
m = 1, . . . , N, the only minimal states are (A, A — 1, . . . , 3, 2, 1) and its cyclic 
shifts (giving A minimal states in all). In this case one has 

M= 

(2!3!---(A-l)!) 2 

An important motivation for much of the previous work on the 2-TASEP 
on Z was its application to the study of shock measures for the 1-TASEP; 
see, for example, [2, 4, 5, 12]. In Section 6, we show how several results 
developed in this context — in particular, concerning the "stationary measure 
as seen from a second-class particle" — can be proved very simply, and often 
strengthened, using the queueing representation. We also give extensions of 
some of these results to the case of the n-TASEP for n > 2. 

In passing we mention another paper [6] in which we discuss (partly in 
survey form) results related to those in this paper, but with a somewhat 
different approach. For example, a proof that product measure is invariant 
for the multiclass TASEP on Z is described much more directly in terms of 
reversibility arguments in infinite volume, in contrast to the combinatorial 
approach in terms of bijections on the state-space of the process on Zjv, used 
in Section 3 of this paper. The aim of [6] is partly to illustrate the key role 
played by reversibility properties, and also to show that the same multiclass 
measures are stationary for a variety of underlying dynamics (including, e.g., 
a multiclass discrete-space version of the Aldous-Diaconis-Hammersley pro- 
cess and certain discrete-time versions of the multiclass TASEP); the method 
of proof is described independently of the dynamics, as far as possible. The 
proof in Sections 3 and 4 of this paper is much more specific to the TASEP, 
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but (perhaps unlike the approach of [6]) adapts quite naturally to the case 
of the process on Zjv- 

2. Main result. 

2.1. n = 2: Angel's construction and the queueing interpretation. We 
first describe Angel's "collapse" process for the case of Z^v- 

Choose pi and p2 with < p\ < p\ + p2 < N. Let A be a subset of 
{1, 2, . . . , iV} of size pi, chosen uniformly at random, and let S be another 
subset, of size P1+P2, chosen uniformly at random and independently of A. 

We will construct a set PCS as follows. T> is initially empty. Consider 
each site i € A in turn (in an arbitrary order). Now we add to the set T> 
the first element in the list {i, i + 1, i + 2, . . .} which is in S and has not yet 
been added to T>. (We "move i to the right," until we find an available space 
for it; the set S is the set of spaces which are available to be used.) Notice 
that a space j E S is "filled" — that is, contained in T> — if and only if there 
is some interval [i,j] which contains at least as many members of A as it 
does members of S (since then, by moving the members of A to the right 
until they encounter a free space, one of them will inevitably be allocated 
to the space j); so indeed, the order in which the sites of A are considered 
does not affect the set T> which is finally obtained by this procedure. Thus 

(2.1) jeV iff jeS and, for some i,\An[i, j]\ >|5n[i, j]\ . 

The set T> constructed in this way (of size p\) gives the locations of first- 
class particles. The set S \ V (of size P2) gives the locations of second-class 
particles. The set S c (of size N — p\ — P2) gives the locations of holes. 

The distribution of the configuration generated in this way is the station- 
ary distribution of the TASEP with 2 types of particle onZjv, with p\ and 
P2 particles of first and second class, respectively, [1]. 

In later notation, we will write, for i £ Z^r, a(i) = 1 if i £ A and 00 oth- 
erwise, s(i) = 1 if i € S and 00 otherwise, and d(i) = 1 if i & T> and 00 

otherwise. Thus a, s and d are configurations from the set u[ = {1, oo} Zjv 
defined earlier. If d is derived from a and s by the process described above 
we will write d = F^ N ' (a, s) . 

The case of Z is analogous. Now we choose Ai and A2 with < Ai < 
Ai + A2 < 1. For each i € Z, let i £ A with probability Ai, independently for 
each i. For each i E Z, let i € S with probability Ai + A2, independently for 
each i and independently of the set A. As in (2.1) we can then define 

(2.2) jeV iff j eS and, for some i < j, \Af) > \S f] [i, 

Equivalently (but rather less formally now that A is infinite), T> can be 
constructed by taking each site in A in turn (in an arbitrary order), and 
moving it to the right until we find an available space in S for it. 
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Again, the configuration of first-class particles (sites of T>), second-class 
particles (sites of S\T>), and holes (sites of S c ) is a sample from the station- 
ary distribution of the TASEP with 2 types of particle on Z, with densities 
Ai and A2 of first- and second-class particles, respectively [1] . 

As above, we will write a(i) = 1 if i G A and 00 otherwise, s(i) = 1 if % G S 
and 00 otherwise, and d(i) = 1 if i G T) and 00 otherwise, where now i G Z. 
Thus a, s and d are configurations from U\ = {l,oo} z . Note that by con- 
struction the random variables a(i),i G Z and s(i),i G Z are all independent. 
If <i is derived from a and s by the process described above we will write 
d = Fx (a, s). 

We now explain the interpretation in terms of a queueing server for the 
case of Z. 

The sites i G Z represent times. The set A represents the set of arrival 
times, when a (single) customer arrives at the queue. The set S is the set of 
service times of the queue. At such times it is possible for a (single) customer 
to depart. The queue thus has arrival rate Ai and service rate fi = X± + A2. 
At each time in S, a customer departs if any customer is present in the 
queue (including one that has just arrived at that same moment). Then T> 
is the set of times at which a customer departs from the queue. S \ T> is the 
set of unused service times. 

We define 



(2.3) Qj 



su P {L4nMI-l«snM|} 

i<3 



Note that since the service rate is larger than the arrival rate, the law of large 
numbers yields that the right-hand side of (2.3) is finite with probability 1. 
From the definition (2.3), and recalling that a(i) = 1 [resp. s(i) = 1] iff i 6 A 
(resp. i 6 S), we can deduce a set of "recurrences," namely that 

(2.4) Qj = + I(a(j) = 1) - I(s(j) = 1)] + 

for each j G Z. Thus Qj may be interpreted as the queue-length after time 
j (since (2.4) has the following interpretation: at each time-step, the queue- 
length increases by 1 if there is an arrival and no service, decreases by 1 if it 
is nonzero and there is a service and no arrival, and otherwise remains the 
same). An equivalent formulation of (2.2) above is that 

(2.5) d(j) = 1 iff s(j) = 1 and either Qj_ x > or a(j) = 1. 

From (2.3), Qj is a function of the variables (a(i), s(i),i < j). Since the 
variables a(i),s(i),i GZ are all independent, (2.4) shows that, conditional 
on the value of Qj, Qj+i is independent of (Q r ,r < j), so that (Qj,j G 
Z) is a discrete-time Markov chain. In fact it is a birth-and-death chain 
with one-step transition probabilities given by i — > i + 1 with probability 
Ai(l — Ai — A2), i — > i — 1 with probability I(i > 0)(1 — Ai)(Ai + A2) and 
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A arrivals 

S services Q # O # 

Fig. 1. Queueing representation of the two-type measure. Arrivals are represented by 
grey particles. Used services (first-class particles) are represented by black particles and 
unused ones (second-class particles) by white particles. 



otherwise i — > i. The stationary distribution of this chain is geometric with 
parameter Ai/(Ai + A2). From (2.3) and the stationarity of a and s, the 
process (Qj,j € Z) is also stationary [and so Qj ~ Geom(Ai/(Ai + A2)) for 
all;]. 

A form of Burke's theorem for the discrete-time M/M/l queue [7] applies 
here. Since the arrival process A is a Bernoulli process with rate Ai, and the 
service process S is a Bernoulli process (independent of the arrivals) with 
rate \i > Ai, the departure process V is itself a Bernoulli process of rate Ai. 

The queueing formalism can also be used to describe the process on Zjv- 
Define 



(2.6) Qj 



su P {\An[i,j]\-\Sn[i,j}\} 



Because \A\ < |«S|, the sup on the RHS of (2.6) is attained by some i which 
satisfies \[i,j] \ < N; that is, i 7^ j + 1 and so the interval is not the 
whole of Ztv- Then, just as before, one can again deduce the relations in 
(2.4) for each j G Z^r, and the equivalence of (2.1) and (2.5) also follows. In 
terms of the collapse process one could give the following description: each 
particle starts at some location i in the process A and is allocated some 
location k in the process S; then Qj is the number of particles for which the 
corresponding i and k satisfy i^k and j 6 [i, k — 1]. 



2.2. Extension to several classes of particle/customer. We now explain 
how to extend this framework to construct the stationary distribution of 
the n-TASEP for any n. In the queueing context, this will correspond to the 
output process of a series of n — 1 queues in tandem with various classes of 
customer; the output process of one queue acts as the input to the next. 

We begin with the case of Z. We first define the behavior of a queueing 
server whose arrival process contains m classes of arrival (and holes) and 
whose departure process contains m + 1 classes of departure (along with 
holes) . 
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Fig. 2. Multi-type representation for 3 types of arrivals. Each arrival has been linked to 
the corresponding service. Note: in this figure and in Figure 3 we have used a "FIFO" 
(first-in first-out) discipline between customer of the same class; when two or more cus- 
tomers present have joint highest priority, the one that departs is the one that arrived ear- 
liest. However, this choice of queueing discipline is entirely arbitrary since for the purposes 
of the TASEP construction, all customers of the same class may be presumed identical. 
For example, the second and the third of the 2nd-class customers in the diagram could 
depart in the opposite order; the departure process d is unchanged. 

Let a = (a(i),i E Z) £ U m be the arrival process, with a(i) €{1,2,..., m}U 
{oo}. The value a(i) = r < oo indicates that a customer of class r arrives at 
time i, while a(i) = oo means that no customer arrives at time i. 

Let s = (s(i),i E Z) E U\ be the service process, with s(i) E {l,oo}. A 
service is available at time i if s(i) = 1, and no service is possible if s(i) = oo. 

We will derive a departure process d = (d(i),i £ Z) £ U m +i from a and s, 
with d(i) £ {1, 2, . . . , m + 1} U {oo}. 

Customers of class 1 have the highest priority, and a customer of class r 
has higher priority than one of class s when r < s. 

If there is a service possible at time i [i.e., s(i) = 1], then the highest 
priority customer of those currently in the queue, including any customer 
who has just arrived at time i itself, departs. (If there is currently more 
than one customer of jointly highest priority, one of these departs.) If the 
departing customer has class k, then d(i) = k. If there is no customer in the 
queue [in particular this implies a(i) = oo], then an "unused service" occurs, 
and we put d(i) = m + 1. (One could imagine that there is a limitless supply 
of customers of class m + 1 available, who get served only when the server 
would otherwise be idle.) If there is no service at time i [i.e., s(i) = oo] then 
put d(i) = oo. No departure occurs; all customers in the queue, along with 
any customer arriving at i, remain in it. Note that d(i) = oo iff s(i) = oo. 

In the situations we need, the arrival and departure processes will be 
stationary and independent, and the long-run intensity of arrivals to the 
queue will be lower than that of services. In this stationary evolution 

for the queue, corresponding to the verbal description just given, can be 
defined in a standard way; see Section 2.3 below. 

We define a function (or "queueing operator") F m :U m x U\ >— > U m+ i, rep- 
resenting the operation of this queue, by F m (a, s) = d. 

Now fix n > 1. Let X be the state-space, {1, oo}Zx{i,...,n} = ^{bV-,"} 
representing configurations of particles on n lines. For x = (x m (i),i £ Z, m £ 
{1, 2, . . . , n}) £ X, we say that x has a particle at site i on line m if x m (i) = 1, 
and that x has a hole at site i on line m if x m (i) = oo. 
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Let X r G (0, 1) for r € {1, 2, . . . , n}, with AH h A n < 1. We consider the 

product distribution on the space X under which all the x m (i) are indepen- 
dent, with ¥(x m (i) = 1) = Ai + h X m . 

Consider a richer state-space V consisting of configurations v = (v m (i), 
iGZ, m £ {1, 2, . . . , n}) with v m (i) 6 {1, 2, ... , m} U {oo}. If v m (i) = oo we 
say that u has a hole at site i on line m, while if v m (i) = r < oo then v has 
a particle of class r at site i on line m. Thus v m , the mth line of v, is a 
configuration in U m ; All the particles on line m have classes in {1,2,..., m}. 

Given a configuration x G A', we will associate a configuration u = G V 
to it. We will have v m (i) = oo iff x m (i) = oo; that is, the particles in v and 
x are in exactly the same places, and deriving v from x corresponds to 
assigning classes to the particles of x. 

First set vi(-) = x±(-). (All the particles on the top line have class 1.) 
Now recursively, define v m+1 (-) = F m (v m (-),x m+ i(-)) for m = 1, 2, . . . ,n - 1, 
where F m is the multiclass queueing operator defined above. Hence we have 
a sequence of n — 1 queues: the mth queue has arrival process v m and service 
process x m+ \ , and outputs the departure process v m+ i (which is then used 
as the arrival process for queue m + 1). 

Note that the intensity of services on line m is Ai + • • • + A m , which is 
increasing in m. Applying Burke's theorem repeatedly, one obtains that for 
r < m < n, the set of particles of class r or lower on line m is a Bernoulli 
process, and in particular has rate (or density) Ai + • • • + A r . The set of 
particles of class exactly r on line m thus has density A r (in the sense that 
each site has probability A r of containing such a particle), but may not be 
a Bernoulli process. 

Theorem 2.1. The distribution ofv m (the mth line of the configuration 
v) is the stationary distribution of the m-type TASEP on Z, with density X r 
of particles of class r. 

There is an analogous result for the case of TL^. First we construct oper- 
ators Fm^ : Um x u[ N ^ i— > U^X analogous to the queueing operators F m 
above. 




Fig. 3. Multi-line representation for 3 types of arrivals. Each arrival has been linked to 
the corresponding service. 
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Suppose we have an "arrival process" (with m types of customer and 
holes) a G Um^ , and a "service process" s G u[ N ^ . We will derive a "departure 
process" (with m + 1 types and holes) d G U^li using an extension of the 
"collapse" described above. 

Let A r = {i G Zjv : a(i) = r} for r <m, and let iS = {i G Zjy : s(i) = 1}. 

Assume that \Ai U • ■ • U A m \ < \S\. We will construct the sets 

T> r = {i: d(i) = r}, r < m + 1. 

First Pi is constructed gradually as follows. T>\ starts empty; take each site 
i G Ai one by one in some arbitrary order, and add to the set T>\ the first 
element of the list + 1, i + 2, . . .} which is in the set <S and which has 
not yet been added to T>\. (As before, the order in which the sites in A\ are 
considered makes no difference to the set T>\ finally obtained.) 

Now suppose that 1 < r < m and that the sets T> s have been constructed 
for s <r. The set T> r is constructed using the same collapse process, based 
on the set A r and the set S \ (T>\ U • • • U XV-i)- Start the set V r empty. 
Take each site i € A r one by one, and add to T> r the first j in the list 
{i, i + 1, i + 2, . . .} such that j E S \ (D\ U • • • U XV- 1) and such that j has 
not yet been added to the set T> r . 

In this way the sets T)\ , . . . , are constructed. Finally let V m+ i = S\ 
(PiU--UD m ). 

Note that \V r \ = \ A r \ for each r <m, and |P m+ i| = |5| — |^4i| |-4. m |- 

We set d(i) = r if i G T> r , for r = 1, . . . , m + 1, and otherwise d(i) = oo. 

The operator is then defined by F^^a,^) = d. 

Fix n> 1. Let be the state-space {1, oo} ZjvX { lv " ,n } of configurations 
of particles with n lines. For x = (x m (i),i G Zat,?tt, G {1,2, . . . ,n}) G X^ N \ 
we say that x has a particle at site i on line m if x m (i) = 1, and that x has 
a hole at site i on line m if x m (i) = oo. 

Consider a richer state space consisting of configurations v = (v m (i), 
i G Z/v, m G {1, 2, . . . , n}) with ^ m (i) G {1, 2, . . . , m} U {oo}. If v m (i) = oo we 
say that v has a hole at site i on line m, while if v m (i) = r < oo then i; 
has a particle of class r at site « on line m. Thus f m , the mth line of v, is a 
configuration in Um^ ; all the particles on line m have classes in {1, 2, ... , m}. 

Given a configuration x G X^ N \ we will associate a configuration v = 
\rW x G to it. We will have ^m(^) = oo iff x m («) = oo; that is, the 

particles in v and x are in exactly the same places, and deriving v from x 
corresponds to assigning classes to the particles of x. 

First set ^i(-) = xi(-). (All the particles on the top line have class 1.) Now 

recursively, define f m +i(-) = F^\v m (-),x m+1 (-)) for m = 1, 2, . . . , n - 1. 

Theorem 2.2. £e£ p r G N /or 1 < r < n, with pi + h p n < N. Let 

1m = Pi + ' ' ' + Pm, /or 1 < m < n. Define X^ to be the subset of X^ 
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consisting of those configurations which have exactly q m particles on line m 
for m = 1,2, . . . ,n. Suppose that x has the uniform distribution on X^\ 
The resulting distribution of the configuration v m (the mth line of the con- 
figuration v) is the stationary distribution of the m-type TASEP on Z^v, 
with p r particles of class r . 



The proofs of Theorems 2.1 and 2.2 are analogous and are done together 
in the next two sections. 

The proof consists of two steps. First we define a Markov process on the 
state-space X (or X^)) (the "multiline process") and show that the prod- 
uct distribution (or the uniform distribution) defined above is a stationary 
distribution for the process. 

Then we show that if x evolves according to the multiline process, then 
v = Vx (or V^x) evolves in such a way that its bottom line (Vx) n realizes 
a TASEP with n types of particle. 

Thus the image of the bottom line of v under the product distribution 
(or uniform distribution) for x is indeed stationary for the n-type TASEP. 

2.3. Further details of multi-class queue construction. In this section we 
give a more formal construction of the queue-length process for the multi- 
class queues described above, and give recurrences analogous to (2.4) for the 
multiclass case. 

First the case of Z,. For k < m, let ^4^j] be the number of customers 
arriving in the interval whose class is less than or equal to k; that 
is, Ajgj = #{r G [i,j]:a(r) < k}. Let S M = #{r e [i,j]:s(r) = 1} be the 
number of services occurring in We now define 

(2-7) Qf = sup[Ajg + . 

The quantity Qj k represents the number of customers in the queue whose 
class is less than or equal to k, just after time j. The definition (2.7) corre- 
sponds to (2.3), applied to the queue in which we consider all customers of 
class less than or equal to k as equivalent, and ignore all other customers. 
Then, analogously to (2.4), 

(2.8) Qf = [Qf\ + I(a(j) <k) — I(s(j) = 1)]. 

The interpretation is as follows: at each time step, the number of customers 
in the queue of class less than or equal to k increases by 1 if such a customer 
arrives and there is no service, decreases by 1 if no such customer arrives 
and there is a service (unless it was already zero), and otherwise stays the 
same. 
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The departure process d = F m (a, s) from the queue is given by: 
for k < to, d(j) < k, 
(2.9) d(j)=m + l, 

d(j) = oo, 

For j G Z, write Qj for the vector (Qj 1 , Qj^, . . . , Qj m )- Note that since 
a(i),s(i),i € Z are all independent, and the processes a and s are stationary, 
we have from (2.7) and (2.8) that the queue-length process Qj,j € Z is 
a stationary Markov chain. Each component Qj k considered in isolation 
corresponds to a simple queue of the sort described in Section 2.1, with 
arrival rate Ai + A2 + • • • + Xk and service rate fx = Ai H h A m+ i . 

In the case of Ztv one defines Auj] and Suj] in the same way and then 

(2.10) Qf = sup[Ag ] -% j] ] + ; 

one again obtains (2.8) and can define d = F^\a,s) by (2.9). For j G Z^r, 
we again write Qj for the vector (Qj , Qj , . . . , Qj m )- It can still be helpful 

to think of Qj k as the "number of customers of type less than or equal to k 
present in the queue just after time j" in order to interpret the recurrence 
(2.8) as the operation of a (priority) queueing server; now this interpretation 
only holds "locally," since the set of "times" is Z^v rather than Z. 

See Figures 1, 2 and 3 for illustrations of the multiclass and multiline 
constructions. 

3. Multiline process. The multiline process with n lines has state space 
X (or XW) = {1, oo}Zx{l,~,»} (or {1, 00 } z ivx{l,...,n} j regp y For m £ {1, . . . , n} 

and i G Z or Ztv, we say that x m (i) is the value at site % on line m. The value 
1 indicates a particle and the value 00 indicates a hole. 

In each of the n lines, the transitions that occur are those of the TASEP 
(with one type of particle): when a bell rings at site i on line to, the configu- 
ration x m (-) moves to xt 1 '^( - )- Bells on the different lines all ring together 
(but not necessarily at the same site). 

We have a process of bells at rate 1 at each site i of the bottom line (line 
n). When a bell rings on line n it generates a bell on the line above (line 
n — 1), which in turns generates one on the line above that (line n — 2) and 
so on up to 1. The rule is as follows. Suppose a bell rings at site i on line 
to. Then if x m (i) = 1 (before any swap), the bell on line to — 1 also rings at 
site i, while if x m (i) = 00, then the bell on line m — \ rings at site i + 1. 

Formally, we can define as follows. For x £ X (or X^ N >) and j G Z (or 
Zjv)) define the sequence b m — 6 m (x, i), to — 0, ... ,77. by b n — i and, for m — 



if s(j) = 1 and either Qj_i > or a(j) < k; 

if s(j) = 1, Qj™ = and a(j) = 00; 
if s(j) = 00. 
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n-l,...,0, 

/g ]\ ^ _(b m+ l, if X rn _|_l(6 m _|_l) = 1, 

m \6 m+ i + l, if x m+ i(b m+ i) = oo. 

Thus b m is the location of the bell which rings on line m, for m = 1, . . . ,n 
(and we have also defined the value bo for later use). So when the bell rings 
at site i on line n, one jumps from the state x to the state y = Y(x, i) defined 
by 

ym{j) = Xrn(j) for j (b rn ~ 1, b m ), 

(3.2) y m {b m - 1) = min{x m (6 m - l),x m (b m )}, 

y m (bm) =max{x m (b m - l),x m (b m )}, 

where b m = b m (x,i) throughout. 

The multiline process is illustrated in Figure 4. 

Theorem 3.1. (\) Let q m e {0,l,...,iV} for me {1,2,. ..,n}. As 6e- 

/ore, /ei Xq be the subset of X^ N } consisting of those configurations which 
have exactly q m particles on line m, for m = 1, 2, . . . , n. Then the uniform 
distribution on is a stationary distribution for the multiline process on 

(\\) Let pi, p2, ■ . ■ , p n € [0,1]. Let u be the product distribution on X under 
which all the {x m (i)} are independent, and ¥(x m (i) = 1) = p m . Then v is a 
stationary distribution for the multiline process on Z. 
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Fig. 4. Multiline process. In this picture we illustrate the effect of two possible bells, at 
i and i' . For i, the bells ring at &4 = 63 = i, 62 = i + 1 and 61 = i + 2, and bo = i + 3. The 
only particle which jumps is on the bottom line, due to the bell at i. For i' the bells ring 
at &4 = 1 and b'- A — b' 2 — b[ = i + 1, and b' — i +2. A particle in the third line jumps due 
to this bell. 
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The following proposition (which is not needed for the main proof) shows 
that, in equilibrium, the multiline process provides a coupling of n station- 
ary TASEP processes (which, by the previous theorem, have independent 
marginals at any fixed time). It will follow very easily from our proof of 
Theorem 3.1; the argument is given at the end of this section. 

Proposition 3.2. Consider the multiline process in one of the station- 
ary distributions described in Theorem 3.1. The marginal law of each line of 
the process is that of a TASEP. 

(In fact, it is possible to go further and show that, in equilibrium, the 
process of bells on any line of the multiline process consists of a collection 
of independent Poisson processes. See [6] for details of this extension in the 
case of Z.) 

The rest of this section is devoted to the proof of Theorem 3.1. We will 
do the proof for (i), the process on Z^v- Part (ii) then follows by letting 
N — ► oo and choosing q[ , . . . ,Qn such that qffl /N — > p m . In that case: 
(a) the uniform measure of part (i) converges to the product measure of 
part (ii); (b) one can construct jointly the processes on Z27V+1 and on Z in 
such a way that the bells at sites N + l,..., 2N, 0, 1, . . . , N of the process on 
Z27V+1 are the same as the bells at sites — N, . . . , — 1, 0, . . . , N of the process 
on Z. Combining (a) and (b), one can obtain a coupling of the sequence of 
finite processes (with initial conditions according to the uniform measure) 
and the infinite process (with initial condition according to the product 
measure) with the following property: for any fixed time t and on any fixed 
finite window [—k, k], the evolution of the processes on Z2N+1 and on Z over 
the time interval [0,t] are identical at all sites in [— k,k] for all large enough 
N, with probability 1. In particular, the expected values of any cylinder 
function of the state at time and at time t for the finite processes converge 
to the corresponding expected values for the infinite process. Hence since 
the uniform measures are stationary for the finite processes, one has also 
that the product measure obtained in the limit is stationary for the infinite 
process. 

So from now on we concentrate on the case of Z^r. For convenience we 
will write X = X^ N) . 

We can describe the process as follows. From any state, the process 
chooses, at rate N, a random i uniformly from {1,2, ... ,N} and jumps 
from its current state x to the state Y(x,i) defined at (3.2). [Note that we 
may have Y(x,i) = x so that the process actually stays in its current state.] 

Consider the function T from X x {1, . . . , N} to itself defined by 



(3.3) 
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where bo is defined at (3.1). 

We will show that T(-, •) is a bijective function. From this it follows that 
the uniform distribution is stationary. To see this, one simply calculates 
the overall rate of jumps into and out of each state y. For convenience we 
include jumps from y to itself which occur under the description above. Let 
fj>(x) = l/|<-f | be the probability of a state x. The rate of jump out of y is 
simply N fj,(y)=N /\X\. For the rate of jump into y, we have 

J2 E m (*) *(y = Y(x, i)) = ^EEE T ( T (x, i) = (y, j)) 

i X I ^ I i X j 

1*1 r 

= JV/I-^I- 

So the rest of the proof is to show that T is indeed a bijective function. We 
will construct its inverse. In fact, this is essentially equivalent to construct- 
ing the time-reversal (in equilibrium) of the process. It turns out that this 
process can be seen as the image on the original process, reflected between 
top and bottom and left and right. 

In the reverse process, the possible transitions on each line are transitions 
of a TASEP in which particles move to the right; that is, (l,oo) may swap 
to (oo, 1). A bell at position j on line m will try to swap positions j — 1 and 
j on line m in this way. Bells now ring at each site on the top line (line 1) 
at rate 1, and each bell on line m generates a bell on line m + 1, as follows. 
If in state y a bell rings at position j on line m, then the bell on line m + 1 
rings at position j if y m (j — 1) = 1 (before any swap) and at position j — 1 
if VmU - 1) = oo. 

Formally, for y G X and j € Z^, define the sequence c m = c m (y,j),m = 
1, . . . , n + 1 by ci = j and, for m = 1, 2, . . . , n, 

•ni if 2/m(Cm 1) — 1) 

m - 1, if y m (c m - 1) = oo. 

Thus when the bell rings at site j on line 1, the bells ring also at sites 
c m (y,j) on lines m = 2, . . . , n, and the process jumps from the state y to the 
state x = X(y,j) defined by 

(3.5) x m (i) = y m (i) for i £ (c m - l,c m ), 

(3.6) x m (c m - 1) = max{y m (c m - l),y m (c m )}, 

(3.7) x m (c m ) = min{y m (c m - 1), y m (c m )}, 



(3.4) 
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Fig. 5. Reverse multiline process. We illustrate the two reverse jumps which correspond 
to the forward jumps in Figure 4. The corresponding reversed bells ring at sites j and j' . 
For j, the bells ring at ci = j, C2 = j — 1, C3 = j — 2 and C4— j — 3; also cs = j — 3. Again, 
only one particle, on the bottom line, jumps due to the bell at j. For j' the bells ring at 
Ci = j' and C2 = C3 = C4 = j' — 1, with c' 5 — j' — 2. Again, a particle in the third line jumps 
due to this bell. 



where c m = c m (y,j) throughout. We have also defined c n+ \. 
This reverse process is illustrated in Figure 5. 
Now define the function T* from X x {1, . . . , N} to itself by 

(3-8) T*(y,j) = (X(y,j),c n+l (y,j)). 

This is the inverse of T that we need: 

Proposition 3.3. IfxeX audi e {1, . . .,N}, thenT*(T(x,i)) = (x,i). 

Proof. Let (y,j) = T(x,i), so that y = Y(x,i) and j = b (x,i). Write 
b m = b m (x,i) and c m = c m (y,j). 

Then the result we need is that X(y,j) = x and c n+ \ = i. 

The configurations x and y are identical except in the following way: for 
each m, if (x m (b m — l),x m (b m )) = (oo, 1) then these values are swapped to 
(l,oo) in y. 

Similarly, the configurations y and X(y,j) are identical except as follows: 
for each m, if (y m (cm — l),ym( c m)) = (1> °°) then these values are swapped 
to (oo,l) in X(y,j). 

So to show that X(y,j) = x, it is enough to show that, for each m, 

(3.9) (x m (b m -l),x m (b m )) = (oo,l) <^> {y m {c m - l),y m (c m )) = (l,oo), 
and that in this case b m = c m . 



(3.14) 



Jm-l 



(3.15) 



Jr-l 
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We first note a few simple implications for later reference. For m = 
l,...,n: 

(3.10) x m (b m ) = 1 => y m (b m - 1) = 1, 

(3.11) £m(frm) = oo =>• y m (b m ) = oo, 

(3.12) (a; m (6 m - l),x m (b m )) = (oo,l) =>• (y m {b m - l),y m (b m )) = (l,oo), 

(3.13) x m (6 m ) = 1 and y m (b m ) = oo =>• x m (6 m - 1) = oo. 

These all follow immediately from the relationship between a; and y observed 
above. 

Now, let C(m) be the property that 

&mj if ^-m(^m) — 1) 

6 m + l, ifx m (6 m ) = oo. 

(The second equality is the definition of 6 m _i.) We will show by induction 
that C(m) holds for m = 1, . . . ,n,n + 1. 

The property C(l) is true by definition, since ci = ci(y,6o) = &o- 

Suppose r G {1, . . . ,n} and C(r) holds, so that 

b r , if x r (b r ) = 1, 

6 r + 1, if x r (b r ) = oo. 

There are two cases to check: 

Case 1. x r (b r ) = 1. Then c r = 6 r [from (3.15)], and also y T (b r — 1) = 1, 
from (3.10). So y r {c T — 1) = 1, and c r+ i = c r = 6 r , using (3.4). Hence C(r + 1) 
holds as desired. 

Case 2. x r (b r ) = oo. Then c r = b r + 1 [from (3.15)], and y r (b r ) = oo 
[from (3.11)]. So y r {c r — 1) = y r (b r ) = oo, and (3.4) gives c r +i = c r — 1 = 6 r , 
giving C(r + 1) as desired. 

Thus by induction C(r) holds for all r = 1, 2, . . . , n + 1. 

In particular putting r = n + 1 we have that c n +i =b n = i. So it remains 
to show that the two sides of (3.9) are equivalent, and that each implies that 

Cm — ^m- 

First suppose that the left-hand side of (3.9) holds, that is, 
(x m (b m - l),x m (b m )) = (oo, 1). 
We have c m = b m by the property (3.14), and then, using (3.12), 
(y m (c m - l),ym(c m )) = (y m (b m - l),2/m(6m)) 
= (l,oo), 

which is the right-hand side of (3.9). In the other direction, suppose that 
the right-hand side of (3.9) holds, that is, 

(y m (c m ~ l),y m (c m )) = (l,oo). 
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Suppose that we also had x m (b m ) = oo. Then we would have y m (b m ) = oo 
from (3.11). Also from (3.14) we would have c m = b m + 1, so that y m (6 m ) = 
ym(c m — l) = 1) which is a contradiction. 

So we must in fact have x m (b m ) = l. Then from (3.14), 

c m — &m as desired. 

We have y m {bm) = 2/m(c m ) = oo, so from (3.13) we have x m (b m — 1) = oo. 
Thus the left-hand side of (3.9) indeed holds as desired. This completes the 
proof of Proposition 3.3 and hence of Theorem 3.1. □ 

Proof of Proposition 3.2. By construction, line 1 of the reverse 
process realizes a (reverse) TASEP. Hence line 1 of the forward process 
realizes a forward TASEP. Now let 1 < m < n and consider the process 
consisting just of lines m, m + 1, . . . , n. This is a multiline process with n — 
m + 1 lines. For exactly the same reason, the top line of this process realizes a 
TASEP. But this is exactly line m of the original process. So (in equilibrium) 
every line realizes a TASEP, as desired. □ 

4. Realization of the multitype TASEP by the multiline process. 

Theorem 4.1 (for both Z and Zjy). If x evolves according to the mul- 
tiline process (with n lines) then v = Vx (or V^ N 'x) evolves in such a way 
that its nth line v n follows a TASEP with n types. 

Since the bells at each site of level n of the multiline process occur at rate 
1 independently, it will be enough to show the following: 

Proposition 4.2. Let m 6 {1,2, . .. ,n}. When a bell rings at site i on 
line m of the multiline process, the effect on v m is also of a bell at site i, 
that is, v m jumps to Vm • 

This is clearly true for m = 1, since v\ =x\. 

For m > 2, recall that a bell at site i on level m creates a bell on level 
m — 1 according to the following rule: if there was a particle at site i on level 
m, then the bell on level m — 1 is also at site i; if not, then the bell on level 
m — 1 is at site i + 1. 

Using this, we can prove the proposition by induction. Under the assump- 
tion that the bell on level m — 1 has the effect of a TASEP bell, we show 
that the swaps on level m — 1 and m also have the effect of a TASEP bell 
on level m. In terms of the arrival, service and departure processes, the in- 
duction step is the following lemma; to perform the induction, we apply it 
with a = v m -i, s = x m and d = v m . 
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Lemma 4.3. (\) Let m>2. Let a be an arrival process with m — 1 types, 
and let s be a service process. Let d = F m (a, s) be the corresponding departure 
process with m types. Define 



+ if s(i) = oo. 

Write d = a( b ~ 1 ' b \ s = s^~ 1 ^ and d = F m (a,s). Then d = d^ 1 ^ . 

(ii) The analogous statement for the process on Zjy rather than Z; the 
operator F m is replaced by . 

Proof. We will word the proof for the case of Z; but an identical argu- 
ment applies to the case of Z^v, using the framework explained at the end 
of Section 2.3. We will check three cases: 

Case 1. s(i) = oo. There is no service at i, and the bell at i leaves the 
process s unchanged: s = s. Also d(i) = oo, so d(* _1 > 1 ) = d. We have b = i + l. 

The order of events in the arrival process at i and i + 1 does not matter; 
since there is no service at i, any arrivals at i or i + 1 have no opportunity 
to depart until i + 1 or later. So d = F m (a, s) = F m (a, s) = F m (a, s) = d = 
c?( l_1 ' 4 ), as required. 

Case 2. s(i — 1) = oo and s(i) = 1. Thus s(i — 1) = s(i) and s(i) = oo. 
Similarly d^~ l,l \i — 1) = d(i) and d^~ l ' % \i) = d(i — 1) = oo. We have b = i. 
Let C be the collection of all the customers who were either present in the 
queue just after time i — 2 or who arrived at time % — 1 or i. The value of 
d(i) is the class of the highest priority customer in the collection C (or is 
equal to m + 1 if C is empty) . 

After the swaps, the higher priority out of the particles arriving at {i — 
1) and i (if any) has moved to place {i — 1) in the arrival process. So the 
particle which previously departed at time i (if any) is now already present 
at time i — Since s(i — 1) = 1, this particle will depart at time i — 1. That 
is, d(i — 1) = d(i), and since s(i) = oo, we have also d(i) = oo = d(i — 1). 
Overall, the collections of particles arriving and departing in [i — are 
unchanged and so Q i , and hence all later values in the queue-length process 
and the departure process, are unchanged. Thus d = d^ % ~ 1 ^ as desired. 

Case 3. s(i — 1) = s(i) = 1. There are services available at both i — 1 and 
i. Consider the collection C as in case 2 above. The values of d(i — 1) and 
d(i) are, in some order, the highest and second highest priorities out of this 
collection C. [The two values may be the same. Also, if this collection has 
fewer than 2 particles, then one or both of d(i — 1) and d{i) will be equal to 
m + 1.] 

The bell has no effect on the service process: s = s. For the arrival process, 
we have b = i, and the swap at {i — 1, i) ensures that d(i — 1) < d(i), so that 
the highest priority customer out of the collection C is available to depart 
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at time i — 1 (since this customer is either already present in the queue 
after time i — 2 or arrives at time i — 1). Thus this customer (which was the 
higher priority of those departing at times % — 1 and i in the process d) will 
now depart at i — 1. The second-highest now departs at i. Again Q i is left 
unchanged. Thus d = d^~ 1,l \ □ 

5. Common denominators and states of minimal weight for the process 
on Zjv 

Theorem 5.1. Consider the n-type process on Zjv, with p m particles of 
type m, for m = 1, 2, . . . ,n, where p± +P2 + • ■ • +p n < N. 

(i) In the stationary distribution, the probability of any state is an in- 
teger multiple of M~ l , where 



(5.1) M ■ 



N\ ( N \ ( N 

Pi) \Pl + P2 J " ' \Pl +P2 H h Pn 



(ii) A state u has probability M~ l (the smallest possible) in the station- 
ary distribution iff the following holds for each j £ Zjy: if u(J) = oo then 
u(j + 1) = oo or u{j + 1) = n, while if u(j) = m < oo then u(j + 1) > m — 1. 

We write for the set of states satisfying the condition of part (ii). 

The proof of the first part is immediate: 

Proof of Theorem 5.1 (i) . M is the number of configurations of the 
multiline process, with (pi + • — \-p m ) particles on line m. By Theorem 2.2, 
the stationary distribution for the n-TASEP is given by the image, under 
the map x — > (V^ N 'x) n , of the uniform distribution for the multiline process. 
Thus the probability of any state is a multiple of M _1 , as required for part 
(i). □ 

Before proving the second part, we need a couple of definitions and lem- 
mas. Recall the "queueing operator" '-Um^ X M[ l— ^m+r ^ ne nex * 

two results concern solutions of the relation d = F^\a,s). Here a £ Um^ 
is an "arrival process" (consisting of m classes of arrival, corresponding 
to values a(j) 6 {1,2, ... ,m}, and holes, corresponding to the value oo); 

is a "service process" (consisting of values 1, where service is pos- 
sible, and values oo, where no service is possible), and d £ is a "depar- 
ture process" (consisting of m classes of departure, corresponding to values 
in {l,2,...,m}, of "unused services," corresponding to the value m + 1, and 
holes or "nonservices," corresponding to the value oo). 

The first lemma shows in effect that s is uniquely determined by d, and 
that there one can always find at least one appropriate a: 
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Lemma 5.2. LetdGli^+i- 

(i) Define H(d) e u[ N) by 

V V " yJ> \l, ifd(j)<oo. 

If (a, s) = d for some a, s, then s = H{d) . 

(ii) Define further G m (d) £ Um^ by 



J oo, if d(j) = m + 1 or d(j) = oo, 

\d(j), ifd(j)<m. 



Then F^\G m (d),H(d))=d. 

Proof. Part (i) follows directly from the definition of Fm\ since if 
d = i* 1 ,™ (a, s), then d(j) = oo iff s(j) = oo. Part (ii) is easy to check. In 
terms of the collapse process, there is a space available directly below the 
starting position of every particle; so each particle remains at the same site 
where it starts. In terms of the queue, there is a service at every time when 
an arrival occurs; thus each arrival departs immediately and the queue is 
always empty. □ 



The following lemma now shows that a = G m (d) gives the unique solution 
henever d G W^_\ 
inherited by G m {d): 



whenever d £ W^_\, and that the property defined in Theorem 5.1(h) is 



Lemma 5.3. If d G W^2i and Fm(a, s) = d, then a = G rn (d) and s = 
H(d). Furthermore, G m (d) eW^fl 

Proof. The fact that s = H(d) is already given by the Lemma 5.2(i), 
so it remains to look at the process a. 

Divide the configuration d into blocks, each of which consists of a string 
of holes followed by a string of particles. Suppose that the block occupies 
the sites io, . . . . . . ,12, with iq < i\ < 22, that d(j) = 00 for j € [tOi*i]j an d 

that d(j) < 00 for j 6 [i\ + l,^]- Note that the definition of W^2i implies 
that d(ii + 1) = m + 1. 

Since there is an "unused service" at i\ + 1, and no services in [io,*i], 
there can be no "arrivals" in [to,*i]; that is, a(j) = 00 for j £ [«0)*i]l and 
the "queue" is empty after time i\. Since there are no more "nonservices" 
in \i\ + 1, 12}, every customer arriving during that time departs immediately. 
Thus for j £ [i\ + 1, 22] , we have a(j) = d(j) if d(j) < m, and a(j) = 00 if 
d(j) = m + 1 (unused service). 
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Hence the block (a(io), . . . , afa)) is identical to the string (d(io), . . . , dfo)), 
except that any value d(j) = m + 1 is replaced by the value a(j) = oo. Putting 
the blocks together, one has exactly that a = G m (d) as asserted. 

Finally, recall that for a configuration in W^_\ the set of "allowable 
pairs" of values for sites j and j + 1 are (oo,oo), (oo,m + 1) or (r,k) for 
some k > r — 1. It's easy to check that replacing all the values m + 1 by 
oo transforms any such allowable pair for W^_\ into an allowable pair for 
Wj£°. Hence G m (d) £ wOP as required. □ 

Proof of Theorem 5.1(h). A state u has probability exactly M~ l iff 
there is a unique x such that v n = u, where v = V^ N 'X. From the definition 
of this means that u m +i = Fk N \v 

mi for each m — 1, 2, . . . , n 1, 

and v\=x\. 

First, we show how to construct one such x directly. Starting from u = v n , 
we can recursively set v m = G m (v m+ i) for m = n — 1, . . . , 1, and x m = H(v m ) 
for each m. Lemma 5.2(h) then gives that indeed v = V^x. 

The configuration v n is a "minimal state" if this is the only possible x. It 
now follows, by applying Lemma 5.3 repeatedly that this is indeed the case 
whenever v n € since then wg can construct v m uniquely from v m+ i 

for each m, and x m uniquely from v m . 

It now remains to show the converse: if u £ VYr/^ , then there are at least 
two x such that {V^ N ^x) n = u. It's enough to show the following: if we are 
given v m+ i £ and are looking for v m such that F^\v m ,x m+ i) = 

w m +i , then either (a) there are two choices for v m or (b) m > 1 and the 
unique choice for v rn gives v m ^ Wm^ ■ This way, starting from u we 
can construct all of v (and hence x) recursively starting from v n = u, and 
be guaranteed to have a nonunique choice at some stage. 

So, choose v m+ i £ W^^. Choose some j such that the pair (v m+ i(j — 1), 

v m+ i(j)) violates the conditions defining w!£2i- 

Case 1. v m+ i(j — 1) = oo and v m+ x(j) < m. No service occurs at j — 1, 
and a customer departs at j. Then v m = G m (v m+ i) satisfies v m (j — 1) = oo 
and v m (j) < m also. But swapping the values v m (j — 1) and v m (j) in the 
arrival process will have no effect on the departure process (since there is 
no service at j — 1). Thus the choice of arrival process v m is nonunique, and 
(a) holds. 

Case 2. u m +i(j) = k < m + 1 and v m+ i(j + l) = r < k — 2. Certainly m > 1. 
Let v m = G m (v m+1 ). Then v m (j) € {k, oo}, and v m (j + 1) <k-2<m-l. 
So v m ^ W%P, and (b) holds. □ 
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6. Independence and comparison properties for the process on Z. Many 
important properties of the 2-TASEP on Z concern the process "as seen from 
a second-class" particle. 

One central result, proved in [4] and [12] may be stated as follows. Let 
u (j)yj € Z be a sample from the stationary distribution (with some densities 
pi and p2 of first- and second-class particles). Then conditioned on u(0) = 
2, the configurations (u(j),j < 0) and (u(j),j > 0) are independent. It is 
sometimes said that the configuration "factorizes" around the position of a 
second-class particle. Put another way, the state 2 is a renewal state for the 
process u(j). 

As observed by Angel [1], there is a simple proof of this fact via the 
collapse process, or, equivalently, via the queueing process. Recall that the 
stationary measure of the 2-TASEP is realized by the departure and service 
process of a single queue. As explained in Section 2, the queueing system is 
Markovian; in particular, given the current queue-length, the future and past 
evolutions of the system are independent (as before, a "time" in the queue 
corresponds to a site of Z for the particle system — the "future" and "past" 
correspond to sites to the right and left respectively). The state u(0) = 2 
corresponds to an unused service at time 0, and therefore implies that the 
queue is empty just after time 0; this information "decouples" the past and 
future as desired. 

For the n-TASEP with re > 3, there is no such single renewal state. How- 
ever, one can find longer "renewal strings." Let r > and let (w(0),w(l), . . . , 
w(r)) be a sequence taking values in {1,2,..., re} satisfying the following 
conditions: 

(i) w(0) = re; 

(ii) w(r) = 2; 

(iii) if n > 4: for each m € {3, .. . ,n — 1}, there exists j m < r such that 
w{jm) = m and such that w(j) < m for all j G {j m + 1, . . . , r}. 

Call such a sequence a renewal string for the re-TASEP. For example, for 
77, = 4, the string (4,1,2,3,1,2) qualifies, with r = 5 and j'3 = 3. One can 
show the following: 

Proposition 6.1. Let w be such a renewal string. Conditional on u(i) = 
w(i) for < i < r, the configurations (u(j),j < 0) and (u(j),j > r) are in- 
dependent. 

Proof. We do not give the proof in full, but the outline is as follows. 
Recall that we can write u = (Vx) n , where x is a sample from the stationary 
distribution of the appropriate multiline process; namely, x is a collection 
of n Bernoulli processes with appropriate rates. We defined a collection of 
re — 1 queueing processes; the mth queue had an arrival process v m (with 
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m — 1 types of customer), service process x m +i, and departure process v m+ \ 
(with m types of customer). Let cj m ^ denote the collection of customers 
present in queue m after time j. Since the queue-length processes form a 
Markov chain as described in Section 2, one has that, given the vector of 
queue-lengths {C^\ ■ ■ ■ , cj™ ^) after time j, the past and future evolution 
of the process v (and in particular its nth line) are independent. Now, the 
key step is to show that if (u(0), . . . ,u(r)) is a renewal string, then Cr is 
empty for all me {1,2, ... ,n— 1}. This can be done by arguments similar to 
those used in the proof of Theorem 5.1; one constructs v m , m = n — 1, . . . , 1 
recursively starting from v n , and uses the information that if a particle of 
type m occurs on line m at site j, then the corresponding (m — l)st queue 
is empty immediately after j. (These arguments are still valid for the case 
of Z, since the proofs of Lemma 5.2 and Lemma 5.3 do not use the cyclic 
structure of Zjv.) From this one obtains that cj is empty for all j m <j<r, 
where js, . . . ,j n -i are defined above and we take j2 = r, j n = 0. □ 

A related property for the 2-TASEP, stated in [2] and proved in [12] and 
[4] may be stated as follows: the set {j > : u(j) = oo} is independent of the 
event {u(0) = 2}. Similarly, by symmetry, the set {j < 0:u(j) = 1} is inde- 
pendent of the event {u(0) = 2}. (The relevant symmetry is the following: 
if one reverses left and right for an n-TASEP and also reverses the order 
of the states {l,2,...,n,oo}, then one again obtains an n-TASEP, with the 
densities appropriately reversed.) 

The following version provides a stronger independence property (and 
also extends to n > 2); it is an immediate consequence of the queueing rep- 
resentation. 

Proposition 6.2. Let u be a configuration distributed according to a 
translation-invariant stationary distribution for the n-TASEP on Z. Then 
the set {j > 0: u(j) = oo} is independent of the entire configuration (u(j) :j < 
0). Similarly, the set {j < 0:u(j) = 1} is independent of the configuration 
«7'),J>0). 

Proof. Write u = (Vx) n as above. The configuration (u(j),j < 0) de- 
pends only on (x m (j),j < 0, m € {1, . . . ,n}), and so is independent of the 
set {j > 0:u(j) = oo}={j > 0:x n (j) = oo} (since the state x consists of n 
independent Bernoulli processes). □ 

In [2, 4] and [12], the densities of particles and holes around a second class 
particle are explored. It is shown that for all j > 



PMj) = l)<P(u(j) = l|u(0) = 2) 
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and similarly (by symmetry), for all j < 

F(u(j) =00) < P(u(i) = oo|it(0) = 2). 

Using the queueing representation, we can extend this to a pathwise stochas- 
tic comparison of the stationary measure with the measure conditioned to 
have a second-class particle at the origin: 

Proposition 6.3. Let u be distributed according to a stationary distri- 
bution for the 2-TASEP on Z ; with positive density of first- and second-class 
particles. Let u' be drawn from the same distribution but conditional on the 
presence of a second-class particle at site 0. Then there is a coupling of u 
and u' such that, with probability 1, one has u(j) < u'(j) for all j > and 
u(j) >u'(j) for all j <0. 

Proof. Let /U < 1 be the rate of the service process S and A < [i the 
rate of the arrival process A, used to construct the stationary distribution 
of the 2-TASEP in Section 2. We have 

Q j = [Q j - 1 +L(jeA)-I(jeS)] + , 

and if we set 

{1, if j G S and either Qj-i > or j G A, 
2, if j G S, Qj-i = and j <£ A, 
00, if j (£ S, 

then u is distributed according to the stationary distribution of the 2-TASEP 
with densities A and fi — A of first- and second-class particles. 

Note that Qo depends only on A(l (— 00, 0] and S n (— 00, 0]; since A and 
S are Bernoulli processes, Qo is independent of An [l,oo) and S H [l,oo). 

Let Qo be drawn from the stationary queue-length distribution for the 
queue (which is geometric with parameter X/fx). 

Under the condition u(0) = 2 we have that the queue-length at time is 
0. So let Q'o = and let Q'j evolve according to 

Q' j = [Q' j - 1 +i(jeA)-l(jeS)] + . 

Set also 

!1, if j G S and either Q'j_\ > or j G A, 
2, if j G <S, Q'j-i = and j £ A, 
00, if j £ S. 

Then (u'(j),j > 0) follows the stationary distribution conditioned on the 
presence of a second-class particle at 0, as desired. Comparing the two sets 
of recurrences, one has Qj > Q'j for all j > 0. The definitions of u and u' 
then give u(j) < u'{j) for all j > as desired. 
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The same argument gives the symmetric statement about the processes 
on j < 0. Since the processes for positive and negative j are independent 
given Qq, the two couplings can in fact be carried out simultaneously as 
claimed. □ 

Note that if T = inf{j > : Qj = 0} then Qj = Q' rj and u(j) = u'(j) for all 
j > T. The distribution of T decays exponentially; thus so does the proba- 
bility that u(j) ^u'(j). Explicit representations for these probabilities can 
be found in [2]. 

We conclude with one more remark concerning the representation of the 
stationary measure of the 2-TASEP by the M/M/l queueing process. We 
assigned first-class particles to the sites of departures and second-class parti- 
cles to the sites of unused services. Using the reversibility of the queue-length 
process, one can obtain that the joint process of arrivals and unused services 
is the time-reversal of the joint process of departures and unused services. 
Thus one could equally assign first-class particles to the sites of arrivals and 
second-class particles to the sites of unused services, to obtain the station- 
ary measure for a 2-TASEP with jumps to the right rather than to the left. 
However, it is less clear how to extend this representation naturally to the 
case n > 2. 

Acknowledgments. We thank the referee for many comments and sug- 
gestions which have substantially improved the presentation of the paper. 
We are grateful for the support of the Fundacao de Amparo a Pesquisa do 
Estado de Sao Paulo (FAPESP), the Conselho Nacional de Desenvolvimento 
Cientffico (CNPq) and the Brazil-France Agreement in Mathematics. 

REFERENCES 

[1] Angel, O. (2006). The stationary measure of a 2-type totally asymmetric exclusion 
process. J. Combin. Theory Ser. A 113 625-635. MR2216456 

[2] Derrida, B., Janowsky, S. A., Lebowitz, J. L. and Speer, E. R. (1993). Ex- 
act solution of the totally asymmetric simple exclusion process: Shock profiles. 
J. Statist. Phys. 73 813-842. MR1251221 

[3] Duchi, E. and Schaeffer, G. (2005). A combinatorial approach to jumping parti- 
cles. J. Combin. Theory Ser. A 110 1-29. MR2128962 

[4] Ferrari, P. A., Fontes, L. R. G. and Kohayakawa, Y. (1994). Invariant measures 
for a two-species asymmetric process. J. Statist. Phys. 76 1153-1177. MR1298099 

[5] Ferrari, P. A., Kipnis, C. and Saada, E. (1991). Microscopic structure of travel- 
ling waves in the asymmetric simple exclusion process. Ann. Probab. 19 226-244. 
MR1085334 

[6] Ferrari, P. A. and Martin, J. B. (2006). Multi-class processes, dual points and 
M/M/l queues. Markov Process. Related Fields 12 175-201. MR2249628 

[7] Hsu, J. and Burke, P. J. (1976). Behavior of tandem buffers with geometric input 
and Markovian output. IEEE Trans. Comm. COM-24 358-361. MR0400453 



28 



P. A. FERRARI AND J. B. MARTIN 



[8] Liggett, T. M. (1976). Coupling the simple exclusion process. Ann. Probab. 4 339- 
356. MR0418291 

[9] Liggett, T. M. (1985). Interacting Particle Systems. Springer, New York. 
MR0776231 

[10] Liggett, T. M. (1999). Stochastic Interacting Systems: Contact, Voter and Exclu- 
sion Processes. Springer, Berlin. MR1717346 

[11] Mallick, K., Mallick, S. and Rajewsky, N. (1999). Exact solution of an exclusion 
process with three classes of particles and vacancies. J. Phys. A 32 8399-8410. 
MR1733708 

[12] Speer, E. R. (1994). The two species asymmetric simple exclusion process. In On 
Three Levels: Micro, Meso and Macroscopic Approaches in Physics (C. M. M. 
Fannes and A. Verbuere, eds.) 91-102. Plenum, New York. 

Instituto de Matematica e Estati'stica Department of Statistics 

Universidade de Sao Paulo University of Oxford 

Caixa Postal 66281 1 South Parks Road 

05311-970 Sao Paulo Oxford OX1 3TG 

Brazil UK 

E-MAIL: pablo@ime.usp.br E-MAIL: martin@stats.ox.ac.uk 

URL: http://www.ime.usp.br/~pablo URL: http://www.stats.ox.ac.uk/~martin 



